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Abstract 



The dispersion coefficients Cq, Cg, and Cio for the interactions between H, 
He, and Li are calculated using variational wave functions in Hylleraas basis 
sets with multiple exponential scale factors. With these highly correlated wave 
functions, significant improvements are made upon previous calculations and 
our results provide definitive values for these coefficients. 
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I. INTRODUCTION 



At large separations R, the interaction potential between two neutral atoms can be 
expressed in terms of inverse powers of R, with the leading term R^^ The nature of 

the long-range interaction can be described by the mutual perturbations of instantaneous 
multipoles of individual atoms. The coefficient of the R~^ term comes from an instantaneous 
dipole-dipole interaction and the coefficient of the R~^ term from an instantaneous dipole- 
quadrupole interaction. 

The precise evaluation of the dispersion coefficients between atoms is computationally 
challenging, because it requires a summation over all intermediate states, including the 
continuum. In actual calculations, it is therefore essential to have an adequate representation 
of the whole spectrum of the Hamiltonian. For atomic systems with more than one electron 
the central problem is the inclusion of electron-electron correlations. 

Recently, significant progress has been made in variational calculations for the 
helium and lithium atoms using double and multiple basis sets in Hylleraas coordinates. 
The nonrelativistic energies for helium have been obtained to better than one part in 10^^ 
for the entire singly-excited spectrum, whereas the nonrelativistic energies for the low-lying 
states of lithium are accurate to a few parts in 10^^ to 10^^ [Q. We have also performed 
a high precision calculation for the lithium — > 2^P oscillator strength |^ which has 
been a subject of controversy for many years. Although there have been many calculations 
for the dispersion coefficients between H, He, and Li (see for example, B,[^-[TT||), the results 
involving Li vary over a considerable range. 



Due to the recent progress of ultracold collisions in both theory and experiment [12 



precise forms of long-range interaction potentials between various atoms become very im- 
portant. The purpose of this paper is to present the results of variational calculations in 
Hylleraas basis sets using multiple nonlinear parameters. The use of our highly correlated 
wave functions will improve upon previous calculations and provide more definitive values 
for the dispersion coefficients. 
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The theory of long-range forces is outhned in Sec. ||. The structures of the basis sets 
for one, two, and three electron atomic systems are presented in Sec. |T|. Our final results 
are tabulated and comparisons are made with the various previous calculations. In the 
appendix, a derivation is given of the dispersion coefficients for the Li(5')-Li(P) system. 



II. FORMULATION 

A. Ground state dispersion coefficients 

In this section, we concentrate on interactions between atoms in their ground states. 
Using second-order perturbation theory, the long-range part of interaction between two 
atoms a and b in their ground states can be expanded in terms of a series of inverse powers 
of the separation R 

V — (^\ 
where the coefficients Ce, Cs, and Cio are 

C6=-Gab(l,l), (2) 
TT 

IS IS 

C8 = T^G,b(l,2) + — G,b(2,l), (3) 
Clo = -Gab(l, 3) + -G'ab(3, 1) + -Gab(2, 2) , (4) 

TT TT TT 

where 

POO 

G^h{L,lh)= al{iuj)a]'^{iuj)duj . (5) 



In d^), af^{iu!) is the dynamic 2^"^ polarizability for atom a at imaginary frequency iu. The 
dynamic polarizability can be expressed in terms of a sum over all intermediate states, 
including the continuum (in atomic units throughout): 

«K-) = Et/^ (6) 



with the 2'-pole oscillator strength /^g being defined by 

/S = ^Ko|(^olE^X(?^)l^")P' (7) 

where EnQ = En — Eq, the sum i runs over all the electrons in the atom, \E'o is the ground 
state wave function, Eq is the corresponding ground state energy, and \E'„ is the nth inter- 
mediate eigenfunction with the associated eigenvalue En- An oscillator strength /^q which 
is independent of magnetic quantum number m is obtained by averaging over the initial 
state orientation degeneracy and summing over the final state degeneracy. It is convenient 
to introduce reduced matrix elements through the Wigner-Eckart theorem ITT 



/ L' I L\ 

{iL'M'\Y.r%a{h)\lLM) = i-lf-'n ]{iL'\\Y.rlYimiL). (8) 

i \-M' m M J i 

With the aid of a sum rule for the 3 — j symbols, the oscillator strength can be written 

where Lq is the total angular momentum for the initial state. 
Using the identity 

du 1 



a, b>0 (10) 



IT Jo (a2 + cu2)(62 + cu2) ab{a + b) 
we can recast Eq. (|^) into the equivalent form 

f{la) r{lh) 

n (J / ^ _ _ \^ JiT-O Jn'o (-,-,\ 

'^abl'-a, 9 pa pb / pa I pb \ ' ^^^^ 

^ nn' -^nO-'^n'ol-^riO + ^^n'o) 

where E^^q = E^ — Eq is the excitation energy for atom i and is always positive for the atoms 
in the ground state. The procedure for evaluating Gab(^a, ^b) is to diagonalize Hamiltonian in 
a basis set and sum over all intermediate states directly according to (|12), and a convergence 
study can be done by increasing the size of basis set progressively. 

The long-range part of the interaction between three ground state atoms is not exactly 
equal to the interaction energies taken in pairs. There is an non-additive term which comes 
from the third order perturbation. The leading terms in the expression for the dispersion 
energy of the three-atom system are |jl| 



T/ Cf Cg'^ Z/abc(3cOsgaCOSgbCOsge + l) 

^ahc ~ 6 6 fi / \S ' V / 

'ab 'be 'ea 1,'ab'be'eaJ 

where 6'a, 6'b, and 6^ are the internal angles of the triangle formed by Tab, ^bc and Tea, and 
i^abe is the triple-dipole constant defined by 
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t'abe = — / a\{iLo)oil{iuj)a'{{iuo)dLO . (13) 
vr JO 



B. Excited state dispersion coefficients 

For two like atoms which are not both in their ground states, the perturbation theory 
for calculating the dispersion coefficients was given by Marinescu and Dalgarno ||10[. They 
worked out all the details for evaluating the dispersion coefficients of alkali-metal dimers in 
different excited states within a one-electron model potential formalism. In this work, we 
examine the important case when one lithium atom is in the 2 "^S ground state and the other 
lithium atom is in the 2 state. A detailed derivation for many-electron systems is given 
at the Appendix. 

The zero-order wave function for the Li(S')-Li(P) system can be written as a symmetrized 
product of two individual atomic wave functions 

= ^[^a(^iMi; r)^b(^2M2; p) + /?^a(^iMi; p)^b(^2M2; r)] , (14) 

where r and p represent all the internal coordinates for the two atoms respectively, Li and L2 
are their total orbital angular momenta. Mi and M2 are the associated magnetic quantum 
numbers, and /5 = ±1 describes the symmetry due to Pauli exclusion principle. Following 
||1CI|| , first-order perturbation theory yields the interaction energy 

y«(L2M2;/?) = -gg±i, (15) 

where 

c^L\i.= (16) 
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The first-order energy correction is proportional to (3. Therefore, for two unlike atoms 
(/5 = 0) in the asymptotic region where overlap between two atoms can be neglected, there is 
no first-order correction to the interaction. For two like atoms, however, there is a first-order 
correction to the interaction energy as long as two atoms are in different angular momentum 
states. For the Li(S')-Li(P) system, the interaction is proportional to R~^. The interaction 
between two ground state atoms is always attractive (see (|l])) but the interaction between 
two like atoms in different angular momentum states is equally likely to be attractive and 
repulsive. 

The leading energy correction obtained from the second-order perturbation theory for 
the Li(S')-Li(P) system is 

^^^^ = -%r, (17) 



where 



i?6 



c^^ = y^^hi — (18) 



with 



n.i = |(vl/,(0;r)||^r,Fi(r,)||x(l;r))| 
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x5:G(l,l,l,A,l,M2)Kvl/b(l;p)||EP.>^i(P.)ll^(A;p))r- (19) 

A j 

In (|18|), the summation should exclude one term which gives rise to E^i^ = E^°\ Note that 
Cg^^ is independent of (3. The values of G are listed in Table |. 

It should be noted that the terms with A = 1 in (|T^, which corresponds to transi- 
tions between even parity states and odd parity state P°, are missing in one-electron 
model potential methods The dominant contribution comes from the radiative tran- 
sition between the lowest doubly excited doublet state ls2p2p^P^, which is stable against 
autoionization [Q, to the singly excited doublet state ls^2p^P°. The contribution of the 



A = 1 term to C^^^ is 0.980 82(5) for Ma = ±1, and 0.392 32(2) for Ma = 0. 
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III. CALCULATIONS AND RESULTS 



For the hydrogen atom, the following Sturmian basis set |15] is used to diagonalize the 
Hamiltonian: 

{r'e-'''-/^L^^'^^\(3r)], (20) 

where L^^'"*"^) (/5r) is the generalized Laguerre polynomial and the parameter (3 is chosen to 
be /3 = 2/(/ + 1). This basis set has proven to be numerically stable as the size of basis set 
is enlarged. 



For the helium atom, the basis set is constructed using Hylleraas coordinates [|16 

{x^3k = r\rir\,e-^^--P'-}, (21) 

and the wave functions are expanded from doubled basis sets. The explicit form for the 
wave function is 

^(ri, r2) = A) + «Sjfc Xijk{oi2, (32)] ± exchange , (22) 

ijk 

and i + j + k < Vt. A complete optimization is then performed with respect to the two 
sets of nonlinear parameters ai, (3i, and 0:2, /52- The screened hydrogenic wave function is 
also included explicitly in the basis set. These techniques yield much improved convergence 
relative to single basis set calculations. 

For the lithium atom, the basis set is also constructed in Hylleraas coordinates 

{cPt,^Mu A, It) = r{' rf r^-^^ r^f ^i^^^-^.r^-p^r.-.^r,^ ^ (23) 

where /xj denotes a sextuple of integer powers ji, j2, ja, ji2, j23, and J31, index t labels 
different sets of nonlinear parameters at, A and 7t. Except for some truncations, all terms 
are included such that 

ji + h + h + J12 + i23 + hi < ^ ■ (24) 
The wave function is expanded from the multiple basis sets 



^(ri,r2,r3) = AY,Y.^i^t^t<Pt,t,t{(^t, Pult) 
t t^t 

X (angular function) (spin function) . (25) 

A similar optimization is also performed with respect to all the nonlinear parameters. 

Table |Tl| contains the values of the static polarizabilities ai(0), a2(0), and ^3(0) for H, 
He, and Li in their ground states. Tables |T| and |I^ present the comparison with selected 
previous calculations for He and Li. Using the Sturmian basis sets containing up to 70 terms 
yields the well-known exact results for the H atom. For He, the largest size of basis set for 
the ground state is 504. For the intermediate states, the largest sizes of basis sets are 728, 
733, and 792 respectively for the P, D, and F symmetries. Table |T| shows that our value 
for ai{0) is in perfect agreement with the best previous results of Bishop and Pipin |^ and 
Jamieson et al. |^ within the first 7 digits. However, our value of ai(0), as well as 02(0) and 
0:3(0), has converged to several more significant figures, as indicated by the extrapolation 
uncertainty in parentheses. For Li, with the fixed size of basis set 919 for the ground state. 
Table contains the convergence studies of ai(0) in both length and velocity forms, as the 
number of terms for the intermediate P symmetry is progressively increased. As a further 
numerical check, we calculated ai(0) for the Li atom by solving an inhomogeneous equation, 
using the Dalgar no- Lewis method |^ as follows: 

«K0) = -^(^(2^5)1 , (26) 

where '^i satisfies the inhomogeneous equation 

(/fo - ^0)^1 + (E r\XU^i) - ^i)*(2 '5) = . (27) 

i 

In the above equation, Hq is the unperturbed Li Hamiltonian, Eq is the Li ground state 
energy, and 

E, = {<^{2^S)\Y.r\YUi^)\'^i:^'S)) . (28) 

i 

Eq. ( P7D is solved variationally by expanding \E^i in terms of a Hylleraas basis set. The 
two procedures are equivalent. Using basis sets up to 2136, ai(0) converges monotonically 
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to 164.109 8. Based on these convergence studies, our final extrapolated value of ai(0) is 
164.111(2). We have also done similar convergence studies for 0:2(0) and 03(0). Our results 
for ai{0) and 0:2(0) are in good accord with the results of Pipin and Bishop ||2^. The model 
potential results of Marinescu et al. agree with the present calculations at the 0.05% and 
0.1% levels for 02(0) and 03(0) respectively. 



Table ^ shows our values of the two-body dispersion coefficients Ce, Cs, and Cio for the 
ground state H, He, and Li atoms. For the H-H case, these coefficients can be calculated 
to arbitrarily high precision. Our value of is in complete agreement with the value of 
Margoliash and Meath Comparisons with the previous calculations for He(l ^S')-He(l ^S) 
and Li(2 ^S)-U{2 ^S) are listed in Tables ^ and [VlT] . 

For the He- He case, our Cq and Cg are in excellent agreement with the results of Bishop 
and Pipin 0, but more precise by about three orders of magnitude. As for Cio, a small 
disagreement of about 1 ppm exists. The value of Jamieson et al. for also agrees with 
our value. 

For the Li-Li case, convergence studies for Cg in both length and velocity forms are listed 
in Table |V]. The agreement between the two forms is satisfactory. From Table |Vlli] , it can 
be seen that the result of Stacey and Dalgarno [^] for Cg is in close agreement with the 
present calculation. The model potential results of Marinescu et al. § for Ce, Cs, and Cio 
differ from our calculations at the 0.1% to 0.3% level. 



Table IX lists the triple-dipole constants i^abc for the combinations between three ground 



state atoms H, He, and Li, together with the previous values of Stacey and Dalgarno 
The overall agreement is about 1%. 

Finally, Tables ^ and ^ list values of C3 and and a comparison with the previous 
calculations for the interaction between the ground state Li(2^S') and the excited Li(2^P). 
C3, which is proportional to the square of the resonant dipole matrix element, has recently 
been calculated to high precision [Q. As for Cg, our values agree with the model potential 
calculations of Marinescu and Dalgarno |]TU[ at about the 0.3% level. 
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APPENDIX A: 

In the appendix, we discuss the dispersion coefficients for the Li(S')-Li(P) system. The 
zero-order wave function for this system can be written as a symmetrized product of two 
individual atomic wave functions 

= ^[^a(i:iMi;r)^b(i:2M2;p)+/?^a(i^iMi;p)^b(i:2M2;r)], (Al) 

where r and p represent all the internal coordinates for the two atoms respectively, Li and L2 
are their total orbital angular momenta. Mi and M2 are the associated magnetic quantum 
numbers, and /5 = ±1 describes the symmetry due to Pauli exclusion principle. Following 
T0| , the interaction potential for two neutral atoms is 



00 00 \r 
VlL 



where 



^ = EE^^ (A2) 

1=1 L=l 



Vi, = (-l)^4vr(/, L)-V2 Y: E Kt,rlpfY,,ir,)Y,.,{p,) . (A3) 

ij fj. 



In the above equation, (/, L, . . .) = {21 + 1)(2L + 1) . . ., and 



J + p / U + /i , 



1/2 



(A4) 



1. The first-order energy 

The first-order energy is given by 
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V^^') = ^^i + ^A2 + /3A3 (A5) 



with 



= (^a(i^iMi; r)^b(i^2M2; p)\V\^^{L,Mi- r)^^{L2M2; p)) , (A6) 

A2 = {^,{L^M,;p)ML2M2;r)\V\^^{L,M^,p)ML2M2;r)), (A7) 

A, = (^a(^iMi; r)^b(^2M2; p)\V\^,{L,M,; p)ML2M2; r)) . (A8) 
Substituting ( |A^ ) into (|A6|) , one has 

IL/i i 

x(*b(i^2M2;p)|^p^^yi-^(p,)|vl/b(L2M2;p)). (A9) 

j 

^From the Wigner-Eckart theorem (H), one has 
(vl>,(LiMi;r)|^rX(f.)l^a(i^iMi;r)) 

i 

= (-l)^^-''M ' ^ ^M(^a(Li;r)||Er^l^(rO||^a(Li;r)). (AlO) 
V-Mi p Ml) i 

For Li = 0, the 3 — j symbol is zero when / > 1. Thus, = 0. Similarly, A2 = 0. For ^3, 

after using the Wigner-Eckart theorem, we obtain 

^3 = E ^^^(^ ^)"'^'^rL(^a(i^i; r)|| E^^^ 

lL^l ^ i 

X (vI/b(L2; p) 1 1 E Pt>^^(P.) I l^a(Li; P))^/ , (All) 

where, for Li = 0, 

( Li I L2\ ( L2 L Li 

V-Mi p M2I \-M2 -p Ml 

= -^^—^^lM^L,L2^f^,-M2 ■ (A12) 

Since [|13 
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(L'||1^||L) = (-1)^'-^(L||11||L'), (A13) 

we finally have 

yW(L.M.;/3) = --^, (A14) 

where 

C,ti, = ^^^47r (^^^^^^^j 1(^,(0; r)|| Er^>^..(f.)ll*b(L2; r))^ ■ (A15) 



2. The second-order energy 

Let the complete set of the system be 

{x{LsM,;r)u;{L,Mf,p)} (A16) 

with the energy eig envalue Eif = + Ef^\ The energy for the unperturbed system is 
= + E2^ . According to the second-order perturbation theory, the second-order 
energy is 

|(^(0)|y|x(i:sM,;rV(LiMi;p))|2 



^1? - 



where T can be written as 



-Y.^^ > (A17) 

IT E^^ - ^(0) ^ ' 



T^Bi + B2 + (5B^ (A18) 



with 



Bi^\Y. E(*a(^i^i;r)*b(^2M2;p)|V^|x(i^.M,;r)a;(LtMi;p))2, (A19) 



^2 = ^ E E(*a(i^i^i;p)*b(i^2M2;r)|V^|x(L,M,;rV(LtMt;p))2, (A20) 
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^3= E E(*a(^i^i;r)*b(^2M2;p)|y|x(i^.M,;rMLtMt;p)) 

LsMs LtMt 

X (^,(LiMi; p)^b(^2M2; r)\V\x{L,M,; r)Lj{LtMf, p)) . (A21) 
After using the Wigner-Eckart theorem, we obtain 

^ ILij. I'L'ijl' LsLt MsMt ^ 

X (*,(Li; r) 1 1 E 1 r)) (*b(i^2; p) 1 1 E pfMpj) I P)) 

i j 

x(*,(Li;r)||Erfyr(?^)||x(i^.;r))(*b(i^2;p)||Epf^^'(P.)ll^(^t;p)) 

/ Li I Ls\f L2 L Lt\/L, I' Ls\/ L2 L' LA 
V-Mi /X mJV-M2 -II Mt)\-M, 11' mJ\-M2 -h' Mtj' 
For Li — 0, the product of four 3 — j symbols becomes 

1 / L2 L LA/L2 L' LA 

Defining Gi by 

L', L,, L„ L2, M2) ^ (-1)^+^' (L, Lo-v^ E ^L-^^Zl' 

/ L2 L LA / L2 L' LA 

(A24) 

V-M2 Mt)\-M2 Ms Mtj 

we have 

^1= E p2L.+L+L^+2 ^1 (^^ L',L,,Lt,L2, M2) I (^a(0; r) 1 1 E rl^'YL. (r,) | |x(L.; r)) |^ 

LL'LsLt i 

X (*b(L2;p)||Ept^4P,)lk(^t;p))(*b(L2;p)||Epf^i'(P.)lk(^*;P))- (A25) 

i i 

Consider the leading term of R~^. The only choice is L^ = 1, L = 1, and L' = 1. If another 
atom is in L2 = 1 state, then Lj = 0, 1, and 2. For this case, 

1 



^i = ;^l(*a(0;r)||Er.yi(?.)llx(l;r))| 



2 



X E 1, 1, A, 1, M2) I (*b(i; p) 1 1 E P.-^i(Pi) I k(A; p)) r ■ (A26) 

A j 

Similarly, for B2 with Li = 0, we have 

13 



E j^2ulw'+2 ^2(1, 1', h,L,,L„M,) I (^,(0; p) 1 1 Y: P-'Yl, (p,) | |a;(L,; p))\' 

ll'LsLt j 

X r) 1 1 E r^ii^i) I IX(^^; r)> (*b(i^2; r) 1 1 E ^f^'(?^) I lx(^.; r)) (A27) 



with 



(47r) 



2 



l{lLt + ij ^^^^ 

I . (A28) 

V-M2 Mi mJ\-M2 Mt M, 
For Lt — 1, I — 1, I' — 1. Thus, for the case of L2 = 1, -B2 becomes 

^2 = ;^|(*a(o;p)||Ep.-^i(p.)IKi;p))r 

x^G2(l,l,l,A,l,M2)|(*b(l;r)||E^^^i(rOllx(A;r))r (A29) 

A i 

Finally, for with Li = 0, we have 
1 

JlLs+Lt+L+l'+2 



^3= E r,L.J^L-,v^2 GsiL,l',L,,L„L,,M,) 

Ll'LsLt 



X (*,(0; r) 1 1 Y: t^^Yls \ r)) (*a(0; p) \ \ E P^Yl, {Pj) \ HLf, p)) 

i j 

X (*b(^2;p)|lEpt^4p,)lk(^t;p))(*b(^2;r)||E^f>^KrOllx(^.;r)), (A30) 

j i 

with 



/ L2 L Lt\/ L2 I' LA 

(A31) 

V-M2 M, MtJ\-M2 Mt Ms) 
The only term which contributes is the one with Ls — 1, Lf — 1, I' — 1, and L = 1. For 
the case of L2 = 1, one obtains 

S3 = -^G'3(1,1,1,1,1,M2) 

X (^^,(0;r)||^r,n(r,)||x(l;r))(^a(0;p)||^p,n(p,)||(^(l;^ 

i j 

X (*b(l;p)|lEP.>^i(P.)ll^(l;P))(*b(l;r)||E^^^i(r^)l|x(l;r)). (A32) 



14 



For the 5* state \l/a(0;r), the parity is +1, and for the P state \l/b(l;r), the parity is — 1. 
Since these two states cannot be connected simultaneously to a third parity eigenstate by a 
dipole operator, B-^ is therefore zero. 

For two like atoms the spectra {x{LM]y)} and {uj{LM\ p)} are identical and Bi and 
B2 can be combined. The final expression for the second-order energy correction to the 
Li(S')-Li(P) system is 



y{2) 




(A33) 



where 



(A34) 



St 



with 



n,t = \{^M^)\\J:r^Y,m\xil;m' 



i 



x$:G(l,l,l,A,l,M2)|(^b(l;p)|lEP.>^i(P.)IKA;p))r 



(A35) 




In ( [A3 51 ), G is defined by 



Gil, 1, 1, A, 1, M2) = 1, 1, A, 1, M2) + G2(l, 1, 1, A, 1, M2) . 



(A36) 



It is easy to see that 



6-1(1, 1, 1, A, 1, M2) = 62(1, 1, 1, A, 1, M2) . 



(A37) 



The algebraic coefficients G(1,1,1,A,1, M2) are listed in Table |. 
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TABLES 



TABLE I. The algebraic coefficients G{1, 1, 1, A, 1, M2). 





A = 




A = 1 


A = 2 


M2 = 

M2 = ±1 


81 " 
81 ' 




81 " 

407^2 
81 " 


176 2 
405 " 

152 2 
405 " 


TABLE IL Values of the static polarizabilities ai(0), 0:2(0), and 0:3(0) for the 


ground state H, 


He, and Li atoms. 










System 


ai(0) 




02(0) 


03(0) 


H 


4.5 




15 


131.25 


He 1.38319217440(5) 




2.445 083101(2) 


10.620 328 6(2) 


Li 


164.111(2) 




1423.266(5) 


39 650.49(8) 


TABLE HI. Comparison of static polarizabilities ai(0), 02(6), and 03(0) for He(l^S'). 


Author (year) 


Reference 


ai(0) 


02(0) 


03(0) 


Luyckx et al. (77) 


0] 


1.379 


2.430 


10.48 


Thakkar (81) 





1.38312 


2.443 44 


10.6144 


Bishop and Pipin (93) 


i 


1.383192 


2.445 083 


10.620 360 


Caffarel et al. (93) 


m 


1.382 7 


2.406 6 


10.36 


Jamieson et al. (95) 


i 


1.383192 






Chen (95) 


@ 


1.383 32 






Chen and Chung (96) 





1.383 27 


2.445 66 


10.625 2 


This work 




1.383192174 40(5) 2.445 083101(2) 


10.620 328 6(2) 
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TABLE IV. Comparison of static polarizabilities ai(0), 02(0), and 0:3(0) for Li(2^S'). 



Author (year) 



Reference 



ai(0) 



"2(0) 



«3(0) 



Maeder and Kutzelnigg (79) 
Muszynska et al. (82) 
Pipin and Woznicki (83) 
Pouchan and Bishop (84) 
Miiller et al. (84) 
Knowles and Meath (86) 
Marouhs and Thakkar (89) 
Pipin and Bishop (92) 
Ponomarenko and Shestakov (93) 
Marinescu et al. (94) 
Wang and Chung (94) 
Merawa et al. (94) 
Kassimi and Thakkar (94) 
Laughhn (95) 
This work 
Experiment (74) 



m 



m 



m 



m 



m 



m 



164.3 

163.8 

163.9 

164(2) 

163.7 

165.8 

164.5 

164.1 

165.2 

164.08 
164.8 

164.2(1) 
163.91 
164.111(2) 
164.0(3.4) 



1383 



1486 
1428 
1423 

1424 

1430 



36 795 



36 495 



39 688 



1423.266(5) 39 650.49(8) 
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TABLE V. Convergence of U{2'^S) ai{0) and U{2'^S)-U{2'^S) Cq in length and velocity forms. 



No. of terms ai(0) (length) 


ai(0) (velocity) 


Cg (length) (velocity) 


56 


164.002 


165.218 


1 389.76 1 409.91 


139 


164.048 


164.201 


1391.21 1393.56 


307 


164.082 


164.131 


1392.37 1393.08 


623 


164.095 


164.107 


1 392.80 1 392.92 


1175 


164.105 


164.108 


1393.17 1393.17 


1846 


164.107 


164.108 


1 393.23 1 393.21 




TABLE VI. \^m1u('s o 


f C\\. C,s. Mild Ci[) for two 


ground stale atoms. 


System 


c& 


Cs 


Cio 


H-H 


6.499 026 705 405 839 313 13 


124.399 083 583 622 343 609 59 3 285.828 414 967 421 697 872 5 


He-He 


1.460 97783768(5) 


14.117857340(5) 


183.691070 5(7) 


Li-Li 


1393.39(16) 


83 425.8(4.2) 


73 721(1) X 10^ 


H-He 


2.821343 915 28(6) 


41.836 376162(8) 


871.540 471(1) 


He-Li 


22.507(1) 


1083.16(5) 


72 602.1(1) 


Li-H 


66.536(5) 


3 279.99(2) 


223 016.6(5) 
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TABLE VII. Comparison of Cq, Cs, and do for the He(l ^5)-He(l ^S) system. 



Author (year) 


Reference 




Cs 


Cio 


Luyckx et al. (77) 





1.458 


14.06 


182.16 


Glover and Weinhold (77) 


[H 


1.459 7(55) 






Margohash and Meath (78) 




1.458 






Bartolotti (80) 


m 


1.463 8 


14.131 


183.47 


Thakkar (81) 


@ 


1.460 82 


14.111 8 


183.600 


Rerat et al. (93) 




1.459 3 


13.883 




Bishop and Pipin (93) 


i 


1.460 977 8 


14.117855 


183.69125 


Jamieson et al. (95) 


[l 


1.460 978 






Chen (95) 





1.461 1 


14.120 


183.74 


Chen and Chung (96) 


Q 


1.46106 


14.120 8 


183.765 


This work 




1.460 977837 68(5) 


14.117857340(5) 


183.691070 5(7) 



TABLE VIII. Comparison of Ce, Cs, and Cio for the Li(2 2S')-Li(2 ^5) system. 



Author (year) 



Reference 



lO-^Ce 



lO-^Cs 



lo-'^c 



10 



Stacey and Dalgarno (68) 
Manakov and Ovsiannikov (77) 
Margohash and Meath (78) 
Maeder and Kutzelnigg (79) 
MiiUer et al. (84) 
Bussery and Aubert-Frecon (85) 
Marinescu et al. (94) 
Merawa et al. (94) 
This work 



HI 

m 



m 



m 
m 



m 



1.391 

1.360 

1.387 

1.389 

1.386 

1.383 

1.388 

1.407 8 

1.393 39(16) 



8.089 

7.578 3 
8.324 
8.43165 
8.342 58(42) 



6.901 

4.816 675 
7.365 

7.3721(1) 
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TABLE IX. Values of the triple-dipole constants fabc for the three ground state atoms H, He, 



and Li. 


kJ J* O IjC-ili 




I -L lllo WUi IS. J 




^6 (.-TtCi. pill ) 


H-H-H 




21 fi42 4fi4 ,^1 n fiS.'i q78 S'^S 1 1 

^ . (JiV7 L»*_>tJ l/ 1 O (JtJO X J- 






He-H-H 




8 1 02 240 874 S('2'l 






Hp-Hp-H 




^ 2fi8 nfi4 8Qfi 1 n 






TTp Hp 




1 A7Q "^^R fiOfi ^ 

1.4: / C7 UtJO UUU O^-LJ 






T 1 N N 
iji-n-n 








27fi 


Li-He-H 




89.830(5) 




89.6 


Li-He-He 




29.824(5) 




29.6 


Ijl-L/l-rl 




6133.5(5) 




6.12 X 10^ 


Li-Li-He 




1917.27(5) 




1.91 X 10^ 


Li-Li-Li 




170 595(6) 




1.69 X 10^ 


TABLE X. 


Values of C3 


and Cq for the interaction between Li(2^S) and Li(2^P). 


Mo 





C3 




Ce 





1 


11.000 226(15) 




2 075.05(5) 





-1 


-11.000 226(15) 




2 075.05(5) 


±1 


1 


-5.500113 3(74) 




1406.08(5) 


±1 


-1 


5.500113 3(74) 




1406.08(5) 
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TABLE XI. Comparison of Cq for the Li(2 25)-Li(2 ^P) system. 



Author (year) 



Reference 



Ce{M2 = 0) 



C6(M2 = ±1) 



Konowalov and Fish (83) 
Vigne-Maeder (84) 
Bussery and Aubert-Frecon (85) 
Marinescu and Dalgarno (95) 
This work 



2 100(50) 
2 025 
1927 
2 066 

2 075.05(5) 



1 750(100) 
1374 
1301 
1401 

1406.08(5) 
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